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ABSTRACT 



The asymmetries of solar p-mode lines provide im- 
portant diagnostic clues as to the nature and location 
of the source of mode excitation. At the same time, 
they can also be responsible for systematic errors in 
the determination of mode frequencies. 

I present a matched wave asymptotic analysis which 
leads to an expression for the spectrum valid for arbi- 
trary stratifications and arbitrary source types. From 
this it is possible to deduce the correct asymmetric 
line profile to use in fitting the spectra, along with 
the physical interpretation of the fitted parameters. 

I apply this procedure to low-degree MDI doppler 
data from which I am able to deduce corrections 
to the mode frequencies. A preliminary comparison 
of the asymmetry parameter with numerical calcula- 
tions yields an estimate of the source depth. 



Key words: solar oscillations; line asymmetry; exci- 
tation. 



1. INTRODUCTION 



The discovery of the asymmetry o f solar p-modc lines 
dates back to Duvall et al. (1993), but it is only with 
recent SOHO and GU1N G data that we can now hope 
to see the asymmetry clearly in individual lines. Con- 
sequently it is only recently that we have begun to 
worry about possible systematic errors in frequency 
determinations which might arise as a result of ig- 
noring the asymmetries. Essentially there are three 
questions in which we are interested :- 1) What is 
the correct line profile to use in fitting asymmetric p- 
mode lines? 2) What can we learn about the source 
of excitation from the asymmetries?, and 3) Why do 
the asymmetries appear oppositely signed in velocity 
and intensity data? In this paper I will concentrate 
on (1) and (2) and comment only briefly on (3). 



There have been numerous models of the asymmetry 
dGabriel 19921 ; |Gabricl 1995| ; |Abrams fe Kumar 1996| ) 
from which it is clear that solar oscillations driven 
near the upper turning point will inevitably have 
asymmetric line-profiles whose properties will be sen- 



sitive to the 
discussion of 



ocation and type of th e source. The 
Rast fc Bogdan (1998)| is particularly 



illuminating because they show that the minima or 
troughs in the spectrum are (at least in some circum- 
stances) eigenfrequencies of a solar model truncated 
from above at the location of the source, and calcu- 
lated with an upper boundary condition determined 
by the source type. The relative location of the peaks 
and troughs is then responsible for the asymmetry of 
the lines. This suggests that, if one could determine 
the locations of these troughs, one could analyse the 
frequencies of these minima using standard helioscis- 
mic techniques. For example, the trough-frequencies 
must satisfy a Duvall Law and a variational princi- 
ple. This conceptual framework provides the basis 
for much of the following discussion. 



2. GENERAL ASYMPTOTIC THEORY 



The general asymptotic theory is described in Rosen- 
thai (1998)| (in what follows, R98). The basic prob- 
lem is the solution of the oscillation equation for a 
plane-parallel stratified layer in the form 
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where the local vertical wavenumber is given by 
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(2) 



and the source is assumed to have the form 



A u ,k{z) = a \K(z )\l 5{z-z )+ai\K(z a )\ ^S'(z-z ). 

(3) 

Note that the solution with a\ — is the Green's 
Function and can therefore be used to generate the 
solution for a truly arbitrary source. For brevity I 
here consider only the low-frequency case in which 
the equation has two turning points and the source 
lies wholly above the upper turning point. The inter- 
mediate and high frequency cases, and the solution 
for an arbitrary source variation with depth, are dis- 
cussed in R98. In the low- frequency case the Fourier 
spectrum in the region above the source is found to 
be 
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where 



Kdz. 



(4) 
(5) 



2 



and 



A = / \K 2 \*dz 



(6) 



and Z t and z t are the lower and upper turning points, 
respectively. 



2.1. Properties of the Peaks and Troughs 



The spectrum diverges at frequencies satisfying 



/ j 

Adz = n 

2 



(7) 



This is a well-known expression for the asymptotic 
frequencies of p modes and, as one would expect, it 
does not depend on the source properties. 



is constant across its depth, and the phase of the 
oscillation is constant across the source, there will 
exist a family of nulls corresponding to frequencies 
at which the amplitude of the solution is identically 
zero in the region above the source. The location of 
these nulls is (obviously) strictly independent of the 
stratification in the region above the source, since the 
amplitude is identically zero there. Thus any changes 
to the atmospheric structure as a result, say, of solar- 
cycle effects, would change the mode frequencies but 
leave the trough frequencies fixed. There therefore 
exists the possibility of a solar-cycle dependence in 
the asymmetry. Any such effect is likely to be ex- 
tremely small as the asymmetry is significant only 
for low-frequency modes whose frequencies vary only 
very slightly over the solar cycle. 



2.1.2. Intensity and Velocity Data 



If the ratio ao + ai/a® — a\ is real, for example if 
either ao or a\ is zero, then the spectrum also has 
zeroes whose location is given by 
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z t 
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If A is large, which is to say that the source lies well 
inside the upper evanescent region, the troughs will 
lie close to the corresponding peaks and the asymme- 
try will be large. The sign of the asymmetry clearly 
depends on the source type. A source containing only 
ao (a monopole source) will have a zero lying imme- 
diately above the mode frequency, and will conse- 
quently have a negative asymmetry, as is seen in the 
MDI and GONG doppler data. For a dipole source 
the sign of asymmery is reversed. 

If A and the source amplitudes are functions only 
of frequency then equations (J?]) and (||) are both in 
the form of Duvall Laws, Fp (w/fc) = n(n + a{u>)/oj) 
with the same Duvall function, Fd, but differing 
phase functions a: 



5a — at 



2 e -2A °o + a i 
7r ao — ai 



(9) 



where a t and a m are the phase functions describ- 
ing the troughs and peaks respectively. Evidently 
the sign of da is positive for a pure monopole source 
and negative for a pure dipole, while its magnitude 
depends quite strongly on the source location. The 
higher the source is placed, the smaller the value of 
da, the closer together the peaks and troughs, and 
the greater the asymmetry. I will therefore use 5a as 
a diagnostic for the source. 



2.1.1. A Solar-Cycle Dependence in the 
Asymmetry? 



Formally the WKB theory predicts that the mode fre- 
quencies depend on the strucure between the turning 
points, while the trough frequencies depend on the 
structure between the lower turning point and the 
source location. However, while the former statement 
is only an approximation, the latter statement can 
actually be made exact. If the phase of the source 



Equation (Q) is the solution for the acoustic wave 
field. Nigam et al (1998), Nigam & Kosovichev 
(1998) and Roxburgh & Vorontsov (1997) all explain 
the reversal of asymmetry between doppler and in- 
tensity measurements as being due to the addition 
of correlated noise, by which they mean that the ac- 
tual measured signal, in at least one of the measure- 
ments, is contaminated by the source. This is equiv- 
alent to adding terms proportional to ao and/or a\ 
to the right hand side of equation (0). It is clearly 
therefore extremely difficult to tell the difference be- 
tween a change of source type and the addition of 
correlated noise. In the data analysis below I will 
essentially assume that the doppler data represents 
the true acoustic wave field and that correlated noise 
is absent in the doppler signal. At present it is not 
clear how that assumption should be validated from 
the data. 



3. Application to Data 

3.1. Asymmetric Line Profiles 

Equation (Ji|) suggests the following form for the 
power in the asymmetric peaks: 

Pi (u) = a 2 + a 3 x 

(7T 7T \ 2 

T (v — vq) + — + iaij 

(10) 

which approximates to 

d r \ _ , (x - a 4 ) 2 + (ai - a 5 ) 2 , , 
P\(y)K,a 2 + a z — — 5 ( n ) 
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where 



x = ~{u - v Q ). 



(12) 



In practice I have found that a 5 , which is non-zero 
only if there is a phase- variation across the source, is 
undcrconstrained by the data, so the final preferred 
form for fitting is 



Pity) = a 2 + a 3 + 
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where the parameters to be fitted are ai, a^, a^, 
the eigenfrequency vq and the full-width at half- 
maximum, r. The asymmetry is evidently due en- 
tirely to a4. These various forms for the asymmetric 
profile are all broad ly similar to that proposed by 
Nigam ct al. (1998)| . However, attention should be 
drawn to the tact that the profiles proposed here are 
based on a general solution to the excitation problem 
valid for an arbitrary stratification and source profile, 
and asymptotic forms for the various parameters (in 
particular a^j can be easily deduced by comparison 
with equation (0), to yield a ready theoretical inter- 
pretation of the fitted values. 



3.2. Data Analysis 

The data analysed consists of 360 day time series of 
MDI doppler-velocity spherical-harmonic coefficients 
for degrees t < 10. (The gap-filled and detrended 
data were provided by Jesper Schou.) The deci- 
sion to look at low-degree modes was based on the 
fact these are less affected by leakage from modes 
of nearby degree. The remaining problem to be ad- 
dressed is therefore the m-leakage from modes of dif- 
fering azimuthal order within the same multiplet. For 
this preliminary analysis, the m-leakage issue was 
sidestepped by looking only at sectoral modes for 
which the effect is negligible. Moreover, the effect of 
m-leakage on the asymmetries of prograde and ret- 
rograde modes should be oppositely directed. Hence 
if the fitted values of «4 agree between prograde and 
retrograde modes this provides an a posteriori justifi- 
cation for neglecting m-leakage. 

The time scries were divided up into ten equal months 
and FFT power spectra calculated for each month. 
The results were then averaged over the ten months. 
Fitted a-coefficients were used to estimate mode fre- 
quencies and all peaks of a particular m/l value in 
a given bin in degree and frequency were averaged 
to find a mean line profile. Figure 1 shows a typ- 
ical example for a low-frequency bin. The different 
modes have very similar line profiles and show a clear 
negative asymmetry. 

The peaks were then averaged to obtain a mean line 
profile, which is then fit to the theoretically derived 
profile (Q3|) . Because each peak is an average of sev- 
eral modes over ten monthly spectra I assume that 
Gaussian statistics are valid and the fitting was there- 
fore carried out using a chi-squared minimisation 
with a weighting proportional to the inverse-square 
of the power. Figure 2 shows the result of the fitting 
for the same case as in Figure 1. Evidently the asym- 
metric profile provides a much better fit to the data 
than the symmetric Lorentzian. 



3.3. Frequency Corrections 

Figure 3 shows the difference in fitted frequencies 
between fits made using symmetric and asymmet- 
ric profiles. For the sectoral modes, the difference 
is small at low frequencies, where the line is so nar- 
row that that the asymmetry makes little difference 
to the fitted frequency. At higher frequencies, the 
lines are sufficiently symmetric that the asymmetry 
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Figure 1. Power spectra of all retrograde sectoral modes 
in the bin 1 <£ < 10 and 1.7 < V < 1.9mHz, aligned using 
fitted a- coefficients. 
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Figure 2. Mean profile of the same modes as in Figure 
1. The dotted curve is a fit of a symmetric Lorentzian 
profile and the dashed curve is a fit of the asymmetric 
profile, Eq. 

coefficient is difficult to determine. The frequency 
error appears to decrease again above 2.5mHz, but 
there is a large scatter. The maximum error is about 
0.1/iHz in the 2 — 2.5mHz range. There is no obvi- 
ous systematic variation in the frequency error as a 
function of degree for these low-degree modes. 



3.4. Asymmetry Parameter 

Figure 4 shows the asymmetry parameter which 
is a measure of the separation of the peak and 
trough. The results for prograde and retrograde sec- 
toral modes are extremely similar (see also Figure 5), 
providing an a posteriori justification for the assump- 
tion that they are unaffected by m-leakage. Once 
again, there is no clear systematic variation with de- 
gree. 
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Figure 3. Frequency difference in the sense (asymmetric 
fit) - (symmetric fit) for sectoral modes in the indicated 
t-bins. Results shown are averages of the results for pro- 
grade and retrograde modes. 



100 J 



2S1S6 



1'6S1S10 



10 




1 1 ^ 

1500 



2000 



2500 
v (/x Hz) 



3000 



3500 



troughs for various source depths are converted 
into phase-function differences by multiplying by an 
empirically-determined Duvall function. The mea- 
sured peak-trough differences are also converted 
to phase-function differences by multiplying by the 
same Duvall function. Since each peak fitted is an av- 
erage over several peaks, error estimates can be made 
by bootstrap resampling of the peaks over which the 
averaging is performed. The results shown in Figure 
5 are broadly consistent with a fairly shallow source, 
located not more than about 250km below the solar 
surface, but evidently the results at this stage are 
somewhat preliminary. 
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Figure 5. Peak-trough separation calculated as a differ- 
ence in phase function for modes in the range 0<£< 10. 
The solid curve is for prograde sectoral modes and the 
dashed curve is for retrograde sectoral modes. Error bars 
are 70% confidence limits. Plotted symbols are numeri- 
cal predictions for various source depths below the pho- 
tosphere :- + -85km, *-200km, o-250km and A-500km. 
Note that 85km corresponds to the location of the sub- 
photospheric peak in the superadiabatic gradient in this 
model. 



Figure 4- Asymmetry parameter (separation between 
troughs and nulls) for the same modes as in Figure 3. 



3.6. Line Widths 



3.5. Comparison with Theory 

Numerical peak and trough frequencies are calculated 
as follows. Eigenmode frequencies for a standard so- 
lar model are determined in the usual way by solving 
the oscillation equations with the appropriate causal 
upper boundary condition. The model is then trun- 
cated at some interior mesh point, corresponding to 
the source location, and the oscillation equations are 
solved with the upper boundary condition Sp = 
applied at that point. The frequencies obtaine d cor- 
respond to the nu lls for the case a\ = 0, which Rast 



fc Bogdan (1998) refer to as a "Sp source" and winch 
Rosenthal (iyy8j calls a monopole source. However 
it should be noted that this is n ot the same usage of 
the ter m "monopole" as that of poldreich & Kumar 
'(1990)1, 



The frequency differences between modes and 



For completeness, Figure 6 shows the fitted line 
widths. The widths are consistent with previously 
known results, particularly the dip in linewidth be- 
tween 2.5mHz and 2.8mHz, providing an extra val- 
idation of the fitting procedure. The inclusion of 
asymmetry has no measurable effect on the fitted 
widths. 



4. Conclusions 



The general asymptotic theory of the excitation of 
solar oscillations demonstrates the relationship be- 
tween the physics of the excitation and the line-shape 
of the mode spectrum. Using this information, I have 
shown how MDI data can be used to determine the 
form and location of the excitation source. The re- 
analysis also provides corrections to the mode fre- 
quencies. The consequences of these corrections for 
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Figure 6. Linewidths for prograde sectoral modes in the 
bin < t < 10. The solid curve is for the asymmetric fit 
and the dashed curve is for the symmetric fit. Error bars 
are 70% confidence limits. 
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more detail in Christensen-Dalgaard et al. (1998) 



In the future, this analysis can be greatly refined by 
the use of more data. In particular, by taking into 
account the m-leakage, it will be possible to use in- 
formation from whole multiplets, not just sectoral 
modes as used here. This will allow much tighter 
limits to be placed on the asymmetry parameter and 
its variation with frequency and degree. Crucially it 
will also then be possible to determine any effect the 
asymmetry may have on the determination of multi- 
plet splittings. 

We also need to consider how to implement the fit- 
ting of these asymmetric profiles into the standard 
analysis packages used to determine helioseismic fre- 
quencies from the data. An important issue is likely 
to be the modelling of the background. A linearly 
sloping background could, in theory, mimic an asym- 
metric line profile, which is why I have here chosen 
a flat background. Standard fitting procedures as- 
sume backgrounds with various power-law forms and 
it will be important to determine if the fitted param- 
eters depend significantly on assum ptions about the 



background (Kumar & Basu 1998). 
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